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I. INTRODUCTION 

Thermo capillarity effects, in which the interfacial properties of two phase systems pro- 
vide the dominant driving forces and control flows and motion, are of interest in a variety of 
areas. [|TJ These include two-phase flows, droplet migration and a variety of other phenom- 
ena in a microgravity environment (see, e.g., Refs. 0|| ), and, for example, convection in 
which the Marangoni effect plays a determining role (see, e.g., Ref. [f| and references cited). 
Similar considerations could apply in transport phenomena mediated by solitons in systems 
with nonuniform temperature. Thus our study may be of interest in charge- density wave 
systems and quasi-one dimensional organic polymers |J and other condensed matter 
systems in which solitons play an important role. @ 

To gain a deeper microscopic understanding into these phenomena, it is of interest to 
consider the effect of a temperature gradient on a two phase interface within a coarse-grained 
description. We approach this question from both analytic and numerical directions. The 
analysis presented here essentially addresses the nature of the free energy of a two-phase 
interface (or kink) in a temperature gradient. Accordingly, the dynamics addressed here are 
restricted to a non-conserved order parameter. 

We find interesting features in the behavior. In particular, a temperature gradient couples 
to the the so-called "translation" mode of the interface which is then only weakly clamped by 
the finite size of the system. A perturbative analysis has a vanishing domain of applicability 
in a large system. The structure of the effective free energy as a function of the collective 
coordinate describing the interface position is sensitive to the boundary conditions, but kinks 
not too near the boundaries are seen to move with constant velocity linear in the gradient. 

The remainder of this paper is organized as follows. In Section [TT] the statics of a non- 
conserved order parameter in a temperature gradient are considered, while in Section |.l [ 
effective Langevin dynamics for the collective coordinate describing the interface position 
are derived. A comparison with a full numerical analysis using relaxational (Model A in the 
lexicon of Ref. M) dynamics is included. Section [IV] is reserved for concluding remarks. 

II. STATICS 

We consider a coarse-grained description of system allowing a transition from a single 
phase to a state of two-phase coexistence with a scalar order parameter, 0. The temperature 
is allowed to vary slowly along one particular direction (chosen to be x) inside the sample, 
which is confined in a box < x < L. We consider configurations translationally invariant 
in directions perpendicular to the gradient. For a configuration <p{x) the free energy per 
unit 'area' is taken to be 



L 

dx 
o 
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In such a model the function r 2 (x) can be taken to specify the local temperature difference 
from a reference (critical) temperature. We take the temperature to be slowly varying and 
parameterize 



r(x) = r (l + -jx\ (2.2) 

This parameterization will prove to be convenient for the analysis to follow. Furthermore, 
for Ar/r <C 1, which will be the case of interest, the temperature varies almost linearly 
with distance, with x = being the hotter and x = L the colder ends of the sample. While 
our explicit calculations and simulations are for the specific model described above, general 
features of the results are not expected to depend on the details of the free-energy functional. 
Equilibrium configurations are extrema of the free energy functional ( |2.1|) , which leads 
to the search for solutions to the following non-linear differential equation: 

d ^ - r 2 {x)<j){x) + \<j) 3 {x) = (2.3) 
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dx 

Notice that with the particular choice ( |2.2| ), <p(x) = ±r(i)/vA is an exact solution of ( |2.3| ) 
corresponding to a particular local equilibrium configuration. It proves convenient to remove 
the local equilibrium variation and to introduce a non-linear change of variables 

4>{x) = ^Vizix)) (2.4) 



and the dimensionless parameter 



Ar , , 



with h <C 1. (Variable changes of this type are discussed in Ref. |J.) This parameter 
measures the strength of the temperature gradient. The dependence of z(x) is determined 
by requiring that the coefficient of the d 2 r]/dz 2 in the resulting differential equation for rj(z) 
be unity. This requirement yields the simple relation 

^P-=r(x). (2.6) 

dx 

We furthermore impose the boundary condition z(0) = thus obtaining the new dimension- 
less variable 

z(x) = r x + -(r x) 2 (2.7) 

measuring the 'distance' from the hot wall. With the assumption that Ar/r <C 1 and hence, 
h, hLro <^C 1, the differential equation for rj(z) in terms of the new dimensionless variable z 
becomes 

rj + 3hr] + r] - rf = (2.8) 

where dots stand for derivatives with respect to z. For h = there are well known "kink" 
solutions to this equation (see e.g., Ref. |]IU|). 



The advantage of parametrizing r(x) as in ( |2.2j ) and of the change of variables ( |2.4| ) 
becomes clear. Whereas in the original differential equation (|2.3|) the x dependence of 



r(x) broke explicit translational invariance, the differential equation for i](z) is manifestly 
translational invariant in the variable z. Furthermore, Eq. ( j2.8|) provides a very appealing 
physical interpretation: i](z) describes the trajectory of a particle moving in "time" z in a 
potential 

V(V) = \f ~ \t (2.9) 

damped by constant friction proportional to h. The present work involves analysis of 
Eq. Q. 

Before searching for solutions to the differential equation, we must specify boundary con- 
ditions. Since we are interested in the behavior of interfaces, relevant boundary conditions 
are those compatible with solutions that have one node (the position of the interface). We 
will focus on two sets of boundary conditions. The first, Type I, corresponds to natural 
boundary conditions: 

Type I: <M U = ; **& |„ 1= 0. (2.10) 

The corresponding conditions on r\{z) at z = and z = I = z[x = L) are rj(z = 0,1) = 
—hr\{z = 0,1). Note we always will consider hLr « 1. 

In the boundary conditions of Type II, the order parameter is fixed at the "local equi- 
librium" values, namely 

Type II: 0(0) = ^M ; <f>{L) = ±^| (2.11) 

corresponding to 7/(0) = =pl ; r)(l) = ±1. 



A. Boundary Conditions I 

1. Perturbative solutions 

For these boundary conditions there is surely a solution in which the order parameter 
remains in one phase, with the \q\ ~ 1. We refer to this as the local equilibrium solution, 
which will be a global minimum of the free energy functional. For small h a perturbative 
solution is easily found with ij(z) = —l + 5(z) and S(z) = e~2 hz (Acosh(Wz) + B sinh(W z)) , 
with A, B chosen to satisfy the boundary conditions and W = J2 + |/i 2 . 

Having established the configuration with the lowest free energy for the case of boundary 
conditions (I), we now concentrate on finding a solution with an interface. This solution 
must necessarily have one node, and for very small h, we expect that its behavior away from 
the node will have \rj\ ~ 1. 

For h = 0, the solution r]o(z) may be found by quadratures in terms of the elliptic sine 



function [II]. For h ^ there are no exact solutions to the differential equation, which in 
its original form (|2.3|) may be recognized (for /« < 1) as a Painleve transcendental without 
an exact solution 0. Although this is a formal property, one of its main consequences is 
that the solution is non-singular in a finite interval (the only possible singularities are at 
infinity), and this suggests that a perturbative approach may be feasible, namely letting 



ri(z) = rjoiz) + h Vl (z) + ■ ■ ■ (2.12) 

This expansion will allow us to study the linear response of the interface profile to the 
temperature gradient. 

We now need to find the first order correction 771. This function obeys the linearized 
equation 

Vi + Vi ~ Sillvi = -3/«7o (2.13) 

This shows that the temperature gradient couples to the "translation mode" of the system, 
5i = rjo, which satisfies the homogeneous equation in Eq. ( |2.13| ). The existence of this mode 
yielding zero eigenvalue for the fluctuation operator (£ = d 2 /dz 2 + l—3r]Q) is a direct result of 
translational invariance (see e.g., Ref. [HJ). Were it not for the fact that the domain is finite, 
such a perturbative analysis could not be carried out. In the present case a solution can 
be constructed using elementary methods. From the solution 5i one can construct another 
linearly independent solution of the homogeneous equation, 62, with unit Wronskian. Some 
additional steps are included in Appendix A. As one might expect, near Z = z — 1/2 ~ 
the solution is of the form 

n(Z) = rio(Z) + hat 6 ^- + ■ ■ • « m (Z + ha x ) + ■■■ , (2.14) 

dZ 

where the constant a.\ (see Appendix A) is determined by the boundary conditions, and 
the terms indicated by dots vanish at Z = 0. Clearly the position of the interface has been 
shifted to Z — —ha,\ (i.e., z = {1/2) — hai); it remains to calculate a.\. A perturbative 
analysis will be possible as long as h(X\ is sufficiently small. 

The explicit evaluation of 7/1 is, in general, complicated by the unwieldly elliptic functions 
and elliptic integrals. However, it is simplified in the physically relevant large volume limit. 
Some background is provided in Appendix B. Since we want only one node of the order 
parameter profile, this requires that the half-period of the elliptic function that corresponds 
to the unperturbed solution flB~T| ) becomes very large in the large volume limit. After some 
tedious but straightforward algebra we find that 

«i « ~e^. (2.15) 

64 

Before proceeding it is useful to consider the perturbative correction in the case of small 
amplitudes. That is, one may linearize Eq. (|2.8|) near rj « 0, discarding the cubic term. Al- 
though in this approximation the differential equation describes an (underdamped) harmonic 
oscillator, the linear response field t]i{Z) can be compared to the exact solution. Requiring 
that there be only one node in the interval requires that / = it (the asymptotic limit of the 
elliptic integral K(m); see Appendix B). The first order correction 7/1 agrees with the first 
order term (in h) in the expansion of the exact solution with boundary conditions of Type 
I. The amplitude of the perturbation is of order hi, and the perturbative analysis is reliable 
whenever W< 1, corresponding physically to the gradient being small on the scale of the 
correlation length. Ultimately secular terms destroy the approximation as hi becomes too 
large. The situation is much worse in the non-linear case to which we now return. 



The result contained in Eq. ( |2.15|) is important: the coefficient that determines the trans- 
lation of the interface becomes exponentially large in the large volume limit, signaling the 
breakdown of perturbation theory. The relevant combination can be seen from Eq. ( |2.14| ) to 



be ha.\. Furthermore, the sign of the coefficient shows that the interface is shifted dramati- 
cally towards the colder end. This suggests that the extremum found does not correspond 
to a local minimum of the free energy functional. However, the result is physically reason- 
able. In the very large volume limit the translation mode is only weakly clamped by the 
boundaries, and since, as we argued above, the temperature gradient couples to this mode, 
the result is a large shift of the interface under the perturbation. 

2. Numerical results 

Since exact analytical solutions are not generally available and perturbation theory leads 
to an exponentially divergent linear response to the temperature gradient, we obtained nu- 
merically the profiles for the order parameter using a shooting method to solve the differential 
equation with Type I (natural) boundary conditions given in Eq. ( J2.10 ). 



Figure 1 shows the rj(z) vs z for h = 10 -3 for a system of size / ~ 40. We clearly see that 
the perturbed interface is established near the cold end of the sample, forming a boundary 
layer of about two correlation lengths. (In dimensionless units, the correlation length is 
£ w y/2). From the numerical standpoint and the interpretation of rj(z) as the trajectory 
of a particle rolling down the potential ( |2.9| ) under constant friction /i<l, the reason that 
the extremum configuration has an interface close to the cold end in the large volume limit 
I ^> £ is clear. 

The initial condition for "shooting", that is the value of r](0) for the boundary condition 
f}(z)\ z=0 = —hrj(0), must be such that 77(0) < —1 ; 77(0) > 0. In this manner, the particle 
first has to climb up the potential hill to 77 = — 1 and reach that point with an extremely 
small velocity, remaining for a long time around that point and slowly falling down on the 
other side of the potential maximum. It then increases its velocity, passes through 77 = 
and climbs up the potential hill towards 77 = +1. Because of friction, the velocity reaches 
zero before reaching the top, and the particle turns back; the integration stops when the 
boundary condition is obeyed again, now with negative velocity. On the other hand, if the 
integration were begun with 77(0) > — 1, for a large volume there would only be solutions 
with many nodes corresponding to a higher free energy. A solution with only one node and 
77(0) > —1 will only appear for small volumes, compatible with the solution found in the 
linearized region near 77 « 0. 

Thus in the large volume limit, the particle must begin with 77(0) < —1 (the "antikink" 
solution will begin with 77(0) > 1) with a small upward velocity given by the boundary 
condition and will remain near the local equilibrium region 77 w — 1 for most of the time, 
making a rather quick transition near the cold end z « I. However unlike in the topological 
"kink" case, the solution with one interface is in the same sector in functional space as the 
local equilibrium solution in that it has the same boundary conditions (of Type I). But 
clearly the local equilibrium solution corresponds to the lowest free energy amongst the 
functions with such boundary conditions. The solutions to the differential equations ( |2.8| ) 
and (|2.8|) are indeed extrema of the free energy functional; the solution with one interface 



cannot be a local maximum because adding "wiggles" in the configuration will increase the 
free energy via the derivative terms. 

If the configuration were a local minimum with a free energy higher than the local equi- 
librium solution, then there must be a local maximum that separates the two solutions. 
However there is no evidence of another solution with the same boundary conditions. This 
reasoning leads us to conjecture that the solution with the interface is most likely a saddle 
point of the free energy functional. To prove this conjecture, we would have to study the 
spectrum of fluctuations around this solution with an interface and identify a particular 
direction in function space for which an eigenvalue is negative. In the present case this is 
an extremely difficult problem complicated by the boundary conditions on the solution. 

In the case of zero temperature gradient and in the large volume limit, because of trans- 
lational invariance, a shift of the interface costs negligible free energy (when the interface is 
far from the boundaries). With the temperature gradient, translational invariance is broken 
and there is a profile which extremizes the free energy. Thus our strategy is to propose 
a good trial "kink-like" function parametrized by the position of the interface z c and to 
compute the free energy as a function of this parameter. This is equivalent to treating the 
position of the interface as a "collective coordinate," which is appropriate in the case of 
kinks and identifies the coordinate, z c , as the translational degree of freedom [P^TB|. 



We have found for a wide range of parameters h and I that the numerical solution to the 
differential equation (|2.8|) is very well described by the interpolating function 



r}(z,z c ) « (1 + — =)tanh 



z - Zc 



(2.16) 



V2 

Here z c determines the position of the interface. The accuracy of this fitting function is 
better than 1% in most of the volume, with slightly larger departures of about 2 — 3% near 
the boundaries of the sample, but extremely accurate near the interface. In terms of r](z; z c ) 
and the variable z, the free energy as a function of z c is given by [up to linear order in h 
consistently with our expansion of the differential equation ( |2.8| )] 




iy- dz{l+3h Ji 



2 



dr}(z;z c )\ 1 2 1 



- nV (2; z c) + - A n 0; Zc) 



dz J 2 ' V ' "' 4 

V 2 (0;z c )) (2.17) 

It is clear from this expression that in terms of the field 77 and the variable z, the breakdown 
of translational invariance is in the metric (and for large enough system, weakly from the 
boundaries) . 

Figure 2 shows F[z c ] vs z c for h = 10~ 3 , I ~ 40 obtained by using fl2.16|) ; these values of 



the parameters are the same as those for Figure 1 (Type I boundary conditions). We see that 
F[z c ] has a maximum at z c = z max « 36 whereas the "shooting" numerical integration gives 
the value of the position of the interface (i.e., the node, 77 = 0) at z c = 35.8, giving confidence 
that the full numerical approach yields a solution corresponding to this maximum. Except 
within a few correlation lengths of the boundaries, we find that F[z c ] varies approximately 
linearly with the interface position z c with a slope 

^fel = ( T 4) 2V2h + 0{h 2 ) + ■■■ = 3ha + 0{h 2 ) + ■■■ (2.18) 

dz c \X I 



where <To = rjj2-\/2/3A is recognized as the interfacial free energy (surface tension) for an 
interface when h = 0. We will compare this to our dynamical simulations described below. 

These and other numerical consistency checks between the trial function approach and 
the full nonlinear solution leads us to conclude that the configuration that extremizes the 
free energy with one interface (node) corresponds to a maximum in the functional direction 
corresponding to translations, and is thus interpreted as a saddle point configuration. This 
saddle already exists for the unperturbed case (h = 0) under the boundary conditions of 
Type I, and the situation for small h represents a smooth deformation. 

If this interpretation is correct, there emerges the question as to the identification of 
the thermodynamically different states separated by this saddle. These states should be 
global minima of the free energy functional, because if one were a global and the other 
a local minimum, there should then be additional solutions of the differential equation 
with different free energies. However, as mentioned above, we find only two: the local 
equilibrium solutions (near rj — ±1) and the interface solution. Thus our conclusion is 
that the saddle point separates the thermodynamically different states corresponding to the 
(nodeless) local equilibrium solutions near t] = ±1, which are degenerate. This interpretation 
will be strengthened by the study of the dynamics in the next section. 

B. Boundary conditions of Type II 

1. Perturbative analysis 

These are "topological" boundary conditions that force the order parameter to have 
at least one node, and clearly there is no equivalent of the local equilibrium configuration 
(in which the order parameter maintains the same sign) that is available with boundary 
conditions I. Hence we expect that the single-node solution of the differential equation 
fl2.8| ) with boundary conditions of Type II is thus an absolute minimum of the free energy 
in the space of functions with these boundary conditions. Before analyzing the solution 
numerically, it proves illuminating to study the linear response as in the previous case. For 
sufficiently large volumes an excellent approximation to the unperturbed solution (h = 0) 
is i] (Z) = tanh[Z/\/2]. With this unperturbed solution we can explicitly construct the 
functions b~i(Z) ; 5 2 (Z) and the response field r\\ satisfying r]i(—l/2) = r] 1 (l/2) = (see 
Appendix A). Again because the boundary conditions are symmetric, the coefficient a 2 
vanishes, and we find (for / 3> 1) 

s/2l 

<*i « ir - (2.19) 

Perturbation theory has a vanishing domain, but an important feature is that, in contrast 
with the previous case, now ct\ is positive. This means that in this case the perturbed 
interface lies closer to the hot end of the sample. 

2. Numerical results 

We have used the same numerical scheme to solve the full nonlinear differential equation 
now beginning with 77(0) = — 1 and "shooting" with an initial derivative such that 
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7](l) = 1. Figure 1 also shows the profile for h = 10~ 3 ; Z = 25 (Type II boundary conditions). 

The profile is again easily understood in terms of the particle rolling down the potential 
hill in the presence of friction: the particle has to begin from 77(0) = —1 but with a fairly 
large derivative because of the friction term. If the initial velocity is small then the particle 
does not make it up the hill to reach rj = 1 Thus the initial velocity is fairly large and the 
particle moves very rapidly initially taking a short time to reach rj = 0; hence the interface 
is very close to z — 0, the hot end. Eventually the particle climbs up the potential hill, 
being slowed down not only by the potential but also by the friction. 

One knows that, because of the coupling to the translation mode, the linear response 
must necessarily diverge in the infinite volume limit. Here, for a large but finite system, the 
position of the perturbed interface is very close (about 2-3 correlation lengths) to the hot 
end, far from the unperturbed value at the middle of the system. The sign of the translation 
is correctly predicted by the linear response calculation as is the case for boundary conditions 
of Type I. However, unlike the case of Type I boundary conditions, this solution, as discussed 
above, corresponds to the lowest free energy compatible with the odd boundary conditions 
of Type II. 

III. DYNAMICS 

A. Langevin equation for the collective coordinate 

Below we will present the results of simulations of relaxational dynamics for the motion 
of a two phase interface in a temperature gradient. Before doing so, we derive from the trial 
function ( |2.16| ) the velocity of the interface. 



Relaxational dynamics for this non-conserved order parameter are specified by assuming 
a Langevin equation description 

d<t>(x, r) 5F 

For our present purposes we neglect a noise term and absorb the characteristic relaxation 
rate into the dimensionless time, r. Assuming that there is a very small distortion of the 
profile as a function of time, that is, that the time evolution corresponds to translations of 
the interface, we propose the parametrization 

4>(x,T) = ^(z-z e (r)) (3.2) 

This parametrization leads at once to an equation for the "collective coordinate" z c (t). For 



£ <C z c <C I and h <C 1 and using the trial function ( |2.16| ) we find 

^ = -3r 2 h + O(h 2 ) + --- (3.3) 

(XT 

The interface is predicted to move with constant speed proportional to the gradient; in the 
language of solitons the coefficient would be identified as the kink (linear) mobility. Details 
are compared directly with the results of simulations below. 
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B. Numerical simulation 

We follow the time evolution of the system using the Langevin dynamics, (|3.1| ). These 
dynamics drive the system to a free energy minimum. Rescaling <ft, x and r to eliminate 
inconsequential parameters, we obtain 

d<j)(x,t) d 2 cj)(x,t) 2 3 

— di — = — dx 2 — ^ ^ ^ X ' ^ ~ ^ ^' '' ^ ' 

where we use the same symbols for the rescaled variables for simplicity. The only param- 
eter remaining is h, which corresponds, as above, to a temperature gradient, such that 
temperature decreases with increasing x if h is positive. 



The equation ( |3.4|) was numerically integrated using a simple Euler discretization on a 
one dimensional lattice of 100 nodes with mesh size Ax = 0.1 and time step At = 0.001. 
We considered separately the two types of boundary conditions discussed above. In Type 
I, the order parameter <ft was required to have zero gradient at the boundaries by imposing 
reflecting boundary conditions. In the second case (Type II), the value of the order parameter 
was fixed at the boundaries, with values corresponding to a different phase at each end of 
the sample. The value chosen was the equilibrium order parameter for an isothermal system 
at the temperature corresponding to that point. (This corresponds to r\ = ±1 in the earlier 
parameterization.) The initial conditions used were of the form 



(x,Q) = (1 + hx) tanh 



X — Xq 



(3.5) 



where xq is the initial position of the kink. This function is very close to the actual values 
that the order parameter takes once it enters the dynamic regime, as long as Xq is not too 
close to the boundaries. 

For the open boundary case (Type I), it was observed that the evolution of the system 
consists of the kink being displaced until it disappears at one of the boundaries. The final 
state is then always a single phase in the whole system, and at that point there is no further 
change in the order parameter. This type of solution has been referred to above as the 'local 
equilibrium' configuration. This is expected, since the one-phase configuration constitutes 
the global minimum of the free energy, and it can be reached with the imposed boundary 
conditions. 

For most initial positions the kink moves towards the higher temperature side. This 
is to be expected since near the high temperature side the correlation length is larger and 
the equilibrium order-parameter gap is smaller. These effects decrease the interfacial free 
energy so that the kink can evolve to a lower free energy state by moving toward the hot 
end. However, it is interesting that if the initial kink position is close enough to the cold 
boundary, the kink disappears at that end. Figures 3(a) and 3(b) show the evolution of the 
kink when it is started from two different positions; in all cases the kink disappears at the 
boundary. This can be understood by looking at the free energy of the system as a function 
of the kink position xo{t), defined as the point where the order parameter vanishes. The free 
energy has a maximum at some value of Xq close to the cold end. The larger the temperature 
gradient h, the closer this maximum is to the cold end. Figure 4 shows this free energy as a 
function of the kink position. It was calculated from equation Q2.1|) applied to the the order 
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parameter configuration as the kink evolves. For a given initial condition, we can obtain 
only a portion of the graph as the interface evolves toward one side. However, for different 
initial conditions, the calculated free energy always falls on the same curve, as expected. 

For the case of fixed boundary conditions, i.e., Type II, the free energy has a minimum 
at some point close to the higher temperature side, and therefore the kink moves towards 
that point from any initial position. The free energy for this type of boundary conditions is 
also shown in Figure 4. 

Finally, also shown in Figure 4 is the free energy calculated for configurations of the 
order parameter given by the trial function ( |3.5p for varying xq. As can be seen, the two 
curves obtained dynamically with differing boundary conditions coincide in the region away 
from the boundaries, which shows that the dynamics of the kink will not be affected by the 
boundary conditions until the kink comes within a few correlation lengths of the boundaries. 
The curve obtained from the local temperature solution (|3.5|) is also very close to the others 
in the region away from the boundaries. 

To compare with the analytical results of equation ( [2.18|) , we computed the value of 



dF[xo]/dxo for small values of h, in the region far from the boundaries, which is where the 
approximate analytical results are expected to hold. The numerical and analytical results 
for the slope of the effective free energy of the interfacial configuration agree within 2%. As 
shown in ( |3.3|) the velocity of the interface is expected to be proportional to the gradient. 
Our numerical simulations yield the coefficient within 2-3% of the approximate analytic 
result. We conclude that as long as the system size is sufficiently large (here ten correlation 
lengths is seen to be large enough), the trial function provides a semiquantitative description 
of the statics and relaxational dynamics. 

IV. CONCLUDING REMARKS 

In this short paper we have considered the statics and non-conserved relaxational dynam- 
ics of a two-phase interface in which the system is subjected to a temperature gradient. For 
the statics one has to extremize the free energy; a non-linear change of variables preserves 
translational invariance and yields a representation in which the interface shape becomes 
equivalent to the time trajectory of a ball rolling in a potential, but slowed by friction 
(Eq. (0)). 

It is seen that the temperature gradient couples to the 'translation' mode of the unper- 
turbed (isothermal) interface, so that perturbation theory can be applied to the introduction 
of the temperature gradient only if the mode is clamped by the finite size of the system. 
Even then the clamping is exceptionally weak for large systems, and the 'linear response' is 
divergent, but correctly predicts which side of the sample will contain the interface. In re- 
alistic terms perturbation theory has vanishing domain since the effective coupling becomes 
~ /iexp(Z), where h is a measure of the gradient and / is a measure of the system size, in 
units of the correlation length. Hence full nonlinear solutions and dynamical simulations 
were carried out numerically. 

The equilibrium configurations are sensitive to the boundary conditions. For Type I 
boundary conditions in which the order parameter derivative vanishes at the walls, the 
global equilibrium is for the system to remain in one phase. A single-kink extremum is 
argued to be a saddle point in the space of functions satisfying these boundary conditions. 
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Simulations using relaxational (Langevin) dynamics reveal that an initial interface travels 
to one of the walls and the interface disappears. An interesting feature is that, generally, 
the interface travels toward the hotter wall at approximately constant velocity; however, if 
the initial interface is established close enough to the colder wall, it is removed there. 

For 'topological' boundary conditions (Type II) in which the order parameter is forced 
to be in two different phases at the ends of the system, the equilibrium configuration has 
the interface near the hotter wall, which is consistent with the behavior of the free energy. 

The structure of the free energy in the space of single kink configurations has been 
analyzed approximately using trial functions which should be accurate as long as the system 
size is sufficiently larger than the thermal correlation length. Good numerical agreement is 
found with purely numerical relaxational dynamics. The trial function along with assumed 
relaxational dynamics yields analytic results for the equation of motion of the interface 
position (collective coordinate), also in good agreement with purely numerical simualtions. 

For an order parameter with conserved (say, Model B ||) dynamics the situation is 
different. A configuration with a single interface cannot change much in response to a 
temperature gradient. However, a kink-antikink pair, representing a slab or bubble of one 
phase in the other can move significantly while respecting the conservation. This topic is 
beyond the present scope and will be explored elsewhere. 

Acknowledgments D.B. thanks H. J. de Vega for illuminating discussions, and he is grate- 
ful to the National Science Foundation for support under grant PHY-93-02534 and INT- 
9016254 (Binational Collaboration with Brazil). F.T. has been supported by CNPQ and 
thanks the Department of Physics and Astronomy of the University of Pittsburgh for its 
hospitality. D.J. and J.L. are grateful to the Microgravity Science and Applications Division 
of NASA for support of the work under grant NAG3-1403. 

APPENDIX A: FORMAL SOLUTION OF THE LINEAR RESPONSE 

The perturbation 771 satisfies Eq. ( J2.13J ). Now, r/ is a function of Z = z — 1/2, so 771 
is also a function of this variable. As noted in the text the function 8\{Z) = t)o(Z) is 
an eigenfunction of the second order fluctuation operator with eigenvalue zero. With this 
solution we can construct another linearly independent with unit Wronskian. Thus we find 
the following solutions of the homogeneous equation: 

h(z) = m(z) 

M-wI'wm (A1) 

The zero mode, Si(Z) is a symmetric function around Z = and vanishes (linearly) 
at Z = ±1/2. This motivates the choice of the lower limit in ( |A1| ), since now 5 2 (Z) is 
antisymmetric around Z = and obviously finite at Z = ±1/2. Finally, the solution to 
Eq. ( ETI3D is 



tji(Z) = 3S 1 (Z) [ Z dZ'S 1 (Z')S 2 (Z') - 35 2 (Z) f dZ'SftZ') + aJ^Z) + a 2 5 2 (Z) (A2) 
Jo Jo 
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with constants a\^ to be determined by the boundary conditions. Expanding the set of Type 
I boundary conditions to linear order in h (and accounting for the explicit /i-dependence of 
the variable Z), one finds that the boundary conditions are symmetric around Z = 0, and 
since b^iZ^) is antisymmetric, the coefficient a 2 vanishes. (This is a bonus of the parametriza- 
tion in terms of Z and the choice of lower limit of the integrals in the functions above.) This 
indicates that the position of the interface is shifted; analyzing the behavior near Z ~ 
yields Eq. (fTT§. 



APPENDIX B: UNPERTURBED KINK 

In a finite box the solution of Eq. ( |2.8|) for h = may be found by quadratures in terms 



of the elliptic sine function [11 



rj (z) = r]iSia[u\m 




u = \z — 

m= A> (m) 

where r] (0) = —rji = —r) (l) gives the values at the end points of the interval. The require- 
ment that there is only one node in the interval is equivalent to requiring that the half-period 
of this solution be identified with I. This requirement in turn determines the value of r\i 
from the relation 



K{m) = l 2 f-^ (B2) 

in which K(m) is the elliptic integral of the first kind |llj and is a quarter period of the 
elliptic function sn. This solution is a function of Z = z — 1/2, and it obeys rjo{Z)\z=~i/2 = 
Vo(Z)\z=i/2 = 0. 
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Figure Captions 

Figure 1 Solid line: r)(z) vs z ; h — 10~ 3 ; I « 40 with boundary conditions (I), show- 
ing the extremum solution with the interface near the cold boundary. Dashed line: 
7](z) vsz ; h — 10~ 3 ; I « 25 with boundary conditions (II), showing the extremum 
(global equilibrium) solution with the interface near the hot boundary. 



Figure 2 F[z c ] vs z c for /i = 10 3 ; I « 40 with the parametrization ( [2.16|) . 



Figure 3 (a) Dynamical evolution of an interfacial structure showing the disappearance of the 
interface at the hot boundary, (b) Evolution of structure with initial interface near 
the cold boundary showing the disappearance of the interface there. 

Figure 4 Effective free energy as a function of the collective coordinate (position of the node in 
the order parameter profile) showing, for comparison, the free energy for the hyperbolic 
tangent trial function. In all cases h = 0.002. 
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